Anisotropy of the optical conductivity of high T c cuprates 
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The optical conductivity along (cr^x) and perpendicular (cr zz ) to the planes is calculated assuming 
strong fen -dependence of the scattering rate and the c-axis hopping parameter. The closed analytical 
expressions for a xx (uj) and a zz (ui) are shown to be integrable at low and high frequencies. A large 
and qualitatively different frequency dependence for both polarizations follows directly from the 
model. The expression for a xx (ui) has an effective scattering rate proportional to frequency, and 
can be easily generalized to provide a simple analytical expression, which may replace the Drude 
formula in the case of non-Fermi liquids. 
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The anomalous transport characteristics n u3 the nor- 
mal state of the high T c superconductorsEra, and the 
corresponding quasi-particle properties in the supercon- 
ducting state have been discussed in a series of recent 
papers emphasizing strong anisotropy of quasiparticles 
as a function of position around the Fermi-surfacecTtl. 
Such anisotropics have been used both in the framework 
of non-pEermi-liquid models and Fermi-liquid-i motivated 
modelsErQ. In the models considered in RefsflEl, anoma- 
lous strong scattering of quasi-particles is assumed in 
the Fermi surface regions closest to the saddle-points, 
while weak scattering is assumed in the zone-diagonal 
directions in fc-space along to (jr, it). Stojkovic and 
Pinesu used a frequency and k-dependent self-energy for 
numerical calculations of the in-plane optical properties 
and magnetotransport of the cuprates. Ioffe and Mil- 
lis (I&M)l3 proposed an analytic expression for the k- 
dependent self energy, containing no explicit frequency 
dependence. The resulting expression for the optical 
conductivity corresponds to a Fermi-surface average of 
Lorentzian response functions with a fc-dependent scat- 
tering rate. Interestingly this expression may be 'correct' 
(at least approximately) irrespective whether the under- 
lying model is Fermi-liquid or not. The strong momen- 
tum depedence of the scattering rate is experimentally 
motivated by Angle Resolved Photo Electron SpectraH 
We will see in this paper, that this model provides an ex- 
cellent description of both the out-of-plane and in-plane 
optical properties of high T c superconductors. We take 
the expression of RefB for the Fermi-surface averaged 
optical conductivity as our starting point 
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In this expression 9 is the angle along the Fermi sur- 
face relative to the (%, tt) direction, n = kp/2ird is the 
density of electrons, m xx = hkp/vF is the planar effec- 
tive mass, and d the distance between planes. For j(0) 
isotropic the classical Drude formula is recovered. (Note 
a factor 1/Ahd difference with Eq. 13 of RefB). I&M 



furthermore postulated the following formula for the k- 
dependent scattering rate, motivated by photoemission 
experiments on a variety of high T c superconductors 



7(ft||) = l/V + rsin 2 (20) 



(2) 



Inserting this expression for the scattering rate in Eq. 
0, substituting £ = e 2ie , carrying out the correspond- 
ing contour integral, and using Cauchi's principle value 
method, we obtain the planar conductivity 

= 7==. 7==i = (3) 



yT — kJry/l + Tt — iur 

where a xx ,o = ne 2 T /m ra . In the approximation of 
RefJa, valid close to the zone-diagonals, the term (uj + 
jr) 1 / 2 i s replaced with the constant (iT) 1 / 2 , resulting 
in non-integrable behavior of YLea xx {uj) at high frequen- 
cies. I&M adopted the model assumption, that 1/r 
should have a quadratic temperature dependence, 1/r = 
T 2 /To, as in a classical Fermi-liquid. (But note, that 
only the zone-diagonal directions have a scattering rate 
which becomes zero at low temperature in their model, 
which is anything but Fermi-liquid behavior). The DC- 
resistivity then has the linear temperature dependence 
p = r 1 / 2 T 1 ^ 2 m xx T / (ne 2 ), observed in optimally doped 
cuprates. 

Let us briefly investigate some of the main analytical 
properties of the above result. First we notice, that 
o-(uj) = a(—oj)* is satisfied. Let us now consider the 
frequency dependent scattering rate, using the definition 
1/t*(lu) = cjReer/Imcr, commonly used for the analysis 
of optical spectra of high T c superconductors and heavy 
Fermion systems 
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(1 + Tt) sjl + turf + y / (l + rr) 2 + (wr) 2 



1 + (lot) 2 + J(1 + Tt) 2 + {ujt) 2 



(4) 

At intermediate frequencies, 1/r <C w <C T, the scat- 
tering rate l/r*(w) = lu. The high frequency limiting 
behavior of these quantities gives 
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lim u xx (uj) = 



m xx lo + i(l/r + r/2) 



(5) 



The low frequency limiting behavior corresponds to a 
Drude conductivity 



Z*(0)ne 2 

o-£)(a;) = limff^w = 



(G) 



with the scattering rate r*(0) and the spectral weight 
Z*(0) defined as 

_ t Tt + 2 _ Tq_ l_ 
T 1 j ~ 2Tr + l ~ 2T 2 + 2r + '"' 



Z*(0) 



f Q °° a D {u))dw ^ 
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where the series expansion is relevant for at low tempera- 
tures, where r diverges. Taken together, we see that a xx 
is integrable, and satisfies 



Rea xx (uj)duj = 

2m xx 



(8) 



We see, that for T < v^oT, which is of the order of 1000 
K, 1/t*(0) is proportional to T 2 , and that the low fre- 
quency spectral weight Z*(0) increases linearly with tem- 
perature. It is amusing, that, apart from a factor two, 
the inverse Hall constant derived in Ref.EJ is precisely this 
spectral weight: 1/Rh = Ba xx /a xy — 2eZ*(Q)n. 
In various papers a linear temperature dependence of the 
scattering rate of the low-frequency Drude tail has been 
reported. However, a fit to Hagen-Rubens behavior pro- 
vides the DC conductivity a xx = Z*(0)ne 2 T*(0)/m xx . 
Usually a temperature independent Drude spectral 
weight is imposed for the purpose of minimizing the 
number of fit parameters. However, the 'Drude scat- 
tering rate' resulting from such a fit, is not a faith- 
ful representation of the actual line- width of o\d(u;) 
if Z*(0) < 1. Instead one can employ the fact, 
that \/t*(oj) — wReffufw) /Ima xx (u>) and Z*(u>) — 
Im[— ne 2 / (4m xx iruxj xx (u))] can be obtained from a careful 
analysis of the DC limiting behavior of the experimen- 
tally measured real and imaginary part of the conductiv- 
ity. This should reveal both the quadratic temperature 
dependence of 1/r and the linear temperature depen- 
dence of the low frequency spectral weight. Although this 
behavior has, to my best knowledge, not been mentioned 
explicitly in the literature, it may be contained in t 
experimental infrared data reported by various groupsl 
This aspect certainly requires further scrutiny, and calls 
for high precision optical experiments in the far infrared 
range. 

In Fig. [I] the evolution of a xx (uj) with temperature is dis- 
played, mimicked by varying the parameter 1/r. This be- 
havior is identical to what has been observed, e«pjerimen- 
tally for the in-plane optical conductivityllSl'Bllil. Note 



also, that the total spectral weight under the optical con- 
ductivity curves is the same for all values of 1/r. 
Let us now investigate the behavior of a zz . Recently 
Xiang and Wheatleyd calculated the c-axis penetration 
depth assuming d-wave pairing, and assuming a model 
for the c-axis transport where momentum parallel to 
the planes, fc||,_is conserved. Based on LDA band the- 
oretical resultsE 2 ] they argued, that for HTSC's with a 
simple (non-body centered) tetragonal structure, t± oc 
(cos(fc x a) — cos(fcj,a)) 2 , leading to a T 5 powerlaw at 
temperature. Following the approach used in RefsJ 
this fc-dependence is evaluated as a function of ku around 
the Fermi surface. ■— ■ 
The c-axis conductivity is, in leading orders of t±£3 



3 2 dm T r2n 



ir 2 H 4 



d0 



t±(0) i 



(9) 



A simple hole, or electron-type surface is circular, and 
has the following k\\ -dependence in leading orders of the 
harmonic expansion 



= h sin 2 (26») 



(10) 



Note, that the expression for a zz now contains a weight- 
ing factor proportional to sin 4 (29), which is sharply 
peaked along (0, ±7r) and (±7r, 0). Again the integral can 
be solved using Cauchi's principle value method, with the 
result 



where 
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iujt)/{Tt), and 
(Jzz.o = e' 2 dm X xt'ol '{itTl a T) . I&M needed a value of KT of 
about 0.15 eV to fit a xx (u>) of YBCOlfJ. In fact, if we use 
Eq.||for a xx (oj), a larger value is required, HT l.2eV to 
produce a good fit. Both values place the c-axis-transport 
well inside the dirty limit region, if we adopt Eq. O as 
an expression representative of the c-axis optical conduc- 
tivity in the normal state of high T c superconductors. 
The analytical formula for 1/r* (a;) according to Eq. [ll] is 
a rather lengthy expression, which is not very illuminat- 
ing. To obtain further insight in the effect of a strongly 
anisotropic scattering rate around the Fermi surface on 
the optical properties, in Fig. || both a xx and a zz and the 
corresponding frequency dependent scattering rates are 
presented, adopting H/t = 12AmeV, and HT — 1.24eV 
as model parameters. Clearly a zz is dominated by the 
transport in the regions close to the saddlepoints, result- 
ing in a frequency dependent scattering rate ranging from 
about r/2 at low frequencies, to T at high frequencies. 
The calculated SBeeteal-shape resembles closely the ex- 
perimental dataB^e-El 

At this point it should be noted however, that the model 
assumption for the fey-dependence of t± is only justified 
for cuprates with the simple tetragonal structure, in par- 
ticular Hgl201 or T11201. Materials like LSCO, T12201, 
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and Bi2212, which have been studied more intensively 
during recent years, have a body centered tetragonal 
structure, for which 

tj_ oc (cos(k x a) — cos(k y a)) 2 cos(k x a/2) cos(k y a/2) (12) 

Hole-type Fermi-surfaces cross not only the zone- 
diagonal directions, but also the lines k x — n and k y = n, 
where extra zeros occur. The integration over 9 in Eq. 
^| thus contains eight zero's, along the (0,±7r), (±7r, 0), 
and (±7r, ±7r) directions. The lowest order harmonic ex- 
pansion of the angular dependence in this case is 



tj_(k\\) = t sin 2 (40) 



(13) 



resulting in the following expression for the c-axis optical 
conductivity 



cr zz (uj) 



(1 + 2ft 2 )(l - 8ft 4 - 8ft 2 ) + 16ft 3 (l + ft 2 ) 3/2 

(14) 



The real part of this expression, and the corresponding 
1/t*(uj) are displayed in Fig. 0, using the same parame- 
ters as before. We observe from this Figure, that for the 
body centered sy terns a zz (uj) is somewhat more peaked. 
This results from the fact that contributions from the 
regions close to the saddle points are suppressed in this 
case. 

Eq. [I] for a xx , suggests a form wich is of general use in 
other materials (not necessarily 2-dimensional) to ana- 
lyze the optical conductivity of non-Fermi liquids with a 
simple analytical formula 



a A (u>) = 



4tt (w + i/Ty- 2a (uj + iT) 2a 



(15) 



If we set r — > oo, and consider only u> <C L,_this corre- 
sponds to the expresion derived by Anderson! 2 ] based on 
the Luttinger liquid model for high T c superconductors. 
The notation for the coefficient a was adopted accord- 
ingly. In Anderson's paper Fermi-liquid response is re- 
stored for a — > 0, as in the present case. Let us briefly 
investigate some of the main analytical properties of a a- 
Again oa(u) — cta(—w)* is satisfied. The high frequency 
limiting behavior of a a 



lim <ta(w) 



1 



4tt (1 - 2a) /t + 2aF - iu 

garantees, that a a is integrable, and satisfies 

I 

Re<JA(u)dui = —u) 



At the low frequency end 

w 2 r (l-2a) r -2a 

hm a A {w) - ^ 1 _ M(1 _ 2a)T + 2a/r] 



(16) 



(17) 



(18) 



which corresponds to a narrow Drude peak with an ef- 
fective carrier life-time r*(0) = (1 - 2a)r + 2a/T. The 
frequency dependent scattering rate 1/t*(lu) crosses over 
to a linear frequency dependence for 1/r <C uj <C T, and 
acquires a constant value (1 — 2a)/r + 2aF for lu — > oo. 
Hence the linear frequency dependence of the scatter- 
ing rate is a robust property of the above phenomcno- 
logical expression for a a- To illustrate this point, in 
Fig. H (Ta(w), the corresponding 1/t*(uj) and 1/Z*(u>) 
are displayed for a number of different values oLql. The 
Luttinger liquid analysis of experimental spectra&Ej gave 
a = 0.15 ± 0.05. The analysis of Ioffe and Millis corre- 
sponds to a = 0.25, which is reasonably close. 
In conclusion, the combination of the k-dependent scat- 
tering rate, and the k-dependence of the inter-plane hop- 
ping parameter known from band theorie, leads to a 
quantitative description of the strongly anisotropic op- 
tical properties found in the high T c cuprates. Both 
the anomalous in-plane optical conductivity and the c- 
axis conductivity follow from the assumption of the k- 
dependent scattering rate, which has been recentely pro- 
posed by Stojkovich, PinesD, Ioffe and MillisEI. Closed 
analytical expressions were obtained for the optical con- 
ductivity, which are really quite simple, and were further 
extended and discussed in a phenomenological framework 
of non- Fermi liquids. The present r esul t .-corresponds 
closely to experimental data of (^(w)!!^^ cia, including 
the high frequency cutoff and the low frequency cross- 
over to Drude behaaripr as can be seen from the satura- 
tion, and of (7^(w)BEJTH 
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FIG. 1. Optical conductivity, a xx , using Eq. 3, with pa- 
rameters fir = 1.2AeV and H/t = 0.124meV (solid curves, 
T = 14 K ), H/t = 1.24meV (dashed curves, T = 44 K 
), H/t = YlAmeV (dotted curves, T = 140 K). The tem- 
peratures correspond to 1/r = T 2 /To, adopting the value 
To = 12meV from Ioffe and Millis. Inset: Frequency depen- 
dent scattering corresponding to the conductivities displayed 
in the main panel. 
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FIG. 2. Optical conductivity using the parameters 
h/r = \2AmeV and KT = 1.24eV. Solid curves o xx , using 
Eq. 3. Dashed: <r zz for the simple tetragonal structure using 
Eq. 11. Dashed-dotted: <r zz for the body centered tetragonal 
structure using Eq. 14. Inset: Frequency dependent scatter- 
ing corresponding to the conductivities displayed in the main 
panel. 
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FIG. 3. Optical conductivity using Eq. 15, with the 
parameters h/r = 1.24meV and ftF = 1.24eV, and 
a — 0.8, 0.4, 0.0, —0.4, and — 0.8 Insets: Frequency dependent 
scattering rate and spectral weight factor. 
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